Duplication and neofunctionalization of a horizontally-transferred xyloglucanase as a facet of the red queen co-evolutionary dynamic.

Oomycetes are heterotrophic protists that share phenotypic similarities with fungi, including the ability to cause plant diseases, but branch in a separate and distant region of the eukaryotic tree of life. It has been suggested that multiple horizontal gene transfers (HGTs) from fungi-to-oomycetes contributed to the evolution of plant-pathogenic traits. These HGTs are predicted to include secreted proteins that degrade plant cell walls. This is a key trait in the pathology of many oomycetes, as the plant cell wall represents a primary barrier to pathogen invasion and a rich source of carbohydrates. Many of the HGT gene families identified have undergone multiple rounds of duplication. Using a combination of phylogenomic analysis and functional assays, we investigate the diversification of a horizontally-transferred xyloglucanase gene family in the model oomycete species Phytophthora sojae. Our analyses detect 11 genes retained in P. sojae among a complex pattern of gene duplications and losses. Using a phenotype assay, based on heterologous expression in yeast, we show that eight of these paralogs have xyloglucanase function, including variants with distinct protein characteristics, such as a long-disordered C-terminal extension that can increase xyloglucanase activity. The functional xyloglucanase variants analysed subtend an ancestral node close to the fungi-oomycetes gene transfer, suggesting the horizontally-transferred gene was a bona fide xyloglucanase. Expression of xyloglucanase paralogs in Nicotiana benthamiana triggers distinct patterns of reactive oxygen species (ROS) generation, demonstrating that enzyme variants differentially stimulate pattern-triggered immunity in plants. Mass spectrometry of detectable enzymatic products demonstrates that some paralogs catalyze production of variant breakdown profiles, suggesting that secretion of multiple xyloglucanase variants increases efficiency of xyloglucan breakdown, as well as potentially diversifying the range of Damage-Associated Molecular Patterns (DAMPs) released during pathogen attack. We suggest that such patterns of protein neofunctionalization, and variant host responses, represent an aspect of the Red Queen host-pathogen co-evolutionary dynamic. Significance Statement: The oomycetes are a diverse group of eukaryotic microbes that include some of the most devastating pathogens of plants. Oomycetes perceive, invade, and colonize plants in similar ways to fungi, in part because they acquired the genes to attack and feed on plants from fungi. These genes are predicted to be useful to oomycete plant pathogens because they have undergone multiple rounds of gene duplication. One key enzyme for attacking plant cell wall structures is called xyloglucanase. Xyloglucanase in the oomycetes has undergone multiple rounds of gene duplication, leading to variants including an enzyme with a C-terminal extension that increases activity. Some xyloglucanase variants trigger unique patterns of reactive oxygen species (ROS) in planta, and generate different profiles of cell wall breakdown products - such outcomes could act to mystify and increase the workload of the plant immune system, allowing successful pathogens to proliferate.

Systematic comparison of unilamellar vesicles reveals that archaeal core lipid membranes are more permeable than bacterial membranes.

One of the deepest branches in the tree of life separates the Archaea from the Bacteria. These prokaryotic groups have distinct cellular systems including fundamentally different phospholipid membrane bilayers. This dichotomy has been termed the lipid divide and possibly bestows different biophysical and biochemical characteristics on each cell type. Classic experiments suggest that bacterial membranes (formed from lipids extracted from Escherichia coli, for example) show permeability to key metabolites comparable to archaeal membranes (formed from lipids extracted from Halobacterium salinarum), yet systematic analyses based on direct measurements of membrane permeability are absent. Here, we develop a new approach for assessing the membrane permeability of approximately 10 µm unilamellar vesicles, consisting of an aqueous medium enclosed by a single lipid bilayer. Comparing the permeability of 18 metabolites demonstrates that diether glycerol-1-phosphate lipids with methyl branches, often the most abundant membrane lipids of sampled archaea, are permeable to a wide range of compounds useful for core metabolic networks, including amino acids, sugars, and nucleobases. Permeability is significantly lower in diester glycerol-3-phosphate lipids without methyl branches, the common building block of bacterial membranes. To identify the membrane characteristics that determine permeability, we use this experimental platform to test a variety of lipid forms bearing a diversity of intermediate characteristics. We found that increased membrane permeability is dependent on both the methyl branches on the lipid tails and the ether bond between the tails and the head group, both of which are present on the archaeal phospholipids. These permeability differences must have had profound effects on the cell physiology and proteome evolution of early prokaryotic forms. To explore this further, we compare the abundance and distribution of transmembrane transporter-encoding protein families present on genomes sampled from across the prokaryotic tree of life. These data demonstrate that archaea tend to have a reduced repertoire of transporter gene families, consistent with increased membrane permeation. These results demonstrate that the lipid divide demarcates a clear difference in permeability function with implications for understanding some of the earliest transitions in cell origins and evolution.

Transporter Proteins as Ecological Assets and Features of Microbial Eukaryotic Pangenomes.

Here we review two connected themes in evolutionary microbiology: (a) the nature of gene repertoire variation within species groups (pangenomes) and (b) the concept of metabolite transporters as accessory proteins capable of providing niche-defining "bolt-on" phenotypes. We discuss the need for improved sampling and understanding of pangenome variation in eukaryotic microbes. We then review the factors that shape the repertoire of accessory genes within pangenomes. As part of this discussion, we outline how gene duplication is a key factor in both eukaryotic pangenome variation and transporter gene family evolution. We go on to outline how, through functional characterization of transporter-encoding genes, in combination with analyses of how transporter genes are gained and lost from accessory genomes, we can reveal much about the niche range, the ecology, and the evolution of virulence of microbes. We advocate for the coordinated systematic study of eukaryotic pangenomes through genome sequencing and the functional analysis of genes found within the accessory gene repertoire.

A Genome Sequence Assembly of the Phototactic and Optogenetic Model Fungus Blastocladiella emersonii Reveals a Diversified Nucleotide-Cyclase Repertoire.

The chytrid fungus Blastocladiella emersonii produces spores with swimming tails (zoospores); these cells can sense and swim toward light. Interest in this species stems from ongoing efforts to develop B. emersonii as a model for understanding the evolution of phototaxis and the molecular cell biology of the associated optogenetic circuits. Here, we report a highly contiguous genome assembly and gene annotation of the B. emersonii American Type Culture Collection 22665 strain. We integrate a PacBio long-read library with an Illumina paired-end genomic sequence survey leading to an assembly of 21 contigs totaling 34.27 Mb. Using these data, we assess the diversity of sensory system encoding genes. These analyses identify a rich complement of G-protein-coupled receptors, ion transporters, and nucleotide cyclases, all of which have been diversified by domain recombination and tandem duplication. In many cases, these domain combinations have led to the fusion of a protein domain to a transmembrane domain, tying a putative signaling function to the cell membrane. This pattern is consistent with the diversification of the B. emersonii sensory-signaling systems, which likely plays a varied role in the complex life cycle of this fungus.

A diversified and segregated mRNA spliced-leader system in the parasitic Perkinsozoa.

Spliced-leader trans-splicing (SLTS) has been described in distantly related eukaryotes and acts to mark mRNAs with a short 5' exon, giving different mRNAs identical 5' sequence-signatures. The function of these systems is obscure. Perkinsozoa encompasses a diversity of parasitic protists that infect bivalves, toxic-tide dinoflagellates, fish and frog tadpoles. Here, we report considerable sequence variation in the SLTS-system across the Perkinsozoa and find that multiple variant SLTS-systems are encoded in parallel in the ecologically important Perkinsozoa parasite Parvilucifera sinerae. These results demonstrate that the transcriptome of P. sinerae is segregated based on the addition of different spliced-leader (SL) exons. This segregation marks different gene categories, suggesting that SL-segregation relates to functional differentiation of the transcriptome. By contrast, both sets of gene categories are present in the single SL-transcript type sampled from Maranthos, implying that the SL-segregation of the Parvilucifera transcriptome is a recent evolutionary innovation. Furthermore, we show that the SLTS-system marks a subsection of the transcriptome with increased mRNA abundance and includes genes that encode the spliceosome system necessary for SLTS-function. Collectively, these data provide a picture of how the SLTS-systems can vary within a major evolutionary group and identify how additional transcriptional-complexity can be achieved through SL-segregation.

A light-sensing system in the common ancestor of the fungi.

Diverse light-sensing organs (i.e., eyes) have evolved across animals. Interestingly, several subcellular analogs have been found in eukaryotic microbes.1 All of these systems have a common "recipe": a light occluding or refractory surface juxtaposed to a membrane-layer enriched in type I rhodopsins.1-4 In the fungi, several lineages have been shown to detect light using a diversity of non-homologous photo-responsive proteins.5-7 However, these systems are not associated with an eyespot-like organelle with one exception found in the zoosporic fungus Blastocladiella emersonii (Be).8Be possesses both elements of this recipe: an eyespot composed of lipid-filled structures (often called the side-body complex [SBC]), co-localized with a membrane enriched with a gene-fusion protein composed of a type I (microbial) rhodopsin and guanylyl cyclase enzyme domain (CyclOp-fusion protein).8,9 Here, we identify homologous pathway components in four Chytridiomycota orders (Chytridiales, Synchytriales, Rhizophydiales, and Monoblepharidiales). To further explore the architecture of the fungal zoospore and its lipid organelles, we reviewed electron microscopy data (e.g., the works of Barr and Hartmann10 and Reichle and Fuller11) and performed fluorescence-microscopy imaging of four CyclOp-carrying zoosporic fungal species, showing the presence of a variety of candidate eyespot-cytoskeletal ultrastructure systems. We then assessed the presence of canonical photoreceptors across the fungi and inferred that the last common fungal ancestor was able to sense light across a range of wavelengths using a variety of systems, including blue-green-light detection. Our data imply, independently of how the fungal tree of life is rooted, that the apparatus for a CyclOp-organelle light perception system was an ancestral feature of the fungi.

Nutrient and salt depletion synergistically boosts glucose metabolism in individual Escherichia coli cells.

The interaction between a cell and its environment shapes fundamental intracellular processes such as cellular metabolism. In most cases growth rate is treated as a proximal metric for understanding the cellular metabolic status. However, changes in growth rate might not reflect metabolic variations in individuals responding to environmental fluctuations. Here we use single-cell microfluidics-microscopy combined with transcriptomics, proteomics and mathematical modelling to quantify the accumulation of glucose within Escherichia coli cells. In contrast to the current consensus, we reveal that environmental conditions which are comparatively unfavourable for growth, where both nutrients and salinity are depleted, increase glucose accumulation rates in individual bacteria and population subsets. We find that these changes in metabolic function are underpinned by variations at the translational and posttranslational level but not at the transcriptional level and are not dictated by changes in cell size. The metabolic response-characteristics identified greatly advance our fundamental understanding of the interactions between bacteria and their environment and have important ramifications when investigating cellular processes where salinity plays an important role.

Emergent RNA-RNA interactions can promote stability in a facultative phototrophic endosymbiosis.

Eukaryote-eukaryote endosymbiosis was responsible for the spread of chloroplast (plastid) organelles. Stability is required for the metabolic and genetic integration that drives the establishment of new organelles, yet the mechanisms that act to stabilize emergent endosymbioses-between two fundamentally selfish biological organisms-are unclear. Theory suggests that enforcement mechanisms, which punish misbehavior, may act to stabilize such interactions by resolving conflict. However, how such mechanisms can emerge in a facultative endosymbiosis has yet to be explored. Here, we propose that endosymbiont-host RNA-RNA interactions, arising from digestion of the endosymbiont population, can result in a cost to host growth for breakdown of the endosymbiosis. Using the model facultative endosymbiosis between Paramecium bursaria and Chlorella spp., we demonstrate that this mechanism is dependent on the host RNA-interference (RNAi) system. We reveal through small RNA (sRNA) sequencing that endosymbiont-derived messenger RNA (mRNA) released upon endosymbiont digestion can be processed by the host RNAi system into 23-nt sRNA. We predict multiple regions of shared sequence identity between endosymbiont and host mRNA, and demonstrate through delivery of synthetic endosymbiont sRNA that exposure to these regions can knock down expression of complementary host genes, resulting in a cost to host growth. This process of host gene knockdown in response to endosymbiont-derived RNA processing by host RNAi factors, which we term "RNAi collisions," represents a mechanism that can promote stability in a facultative eukaryote-eukaryote endosymbiosis. Specifically, by imposing a cost for breakdown of the endosymbiosis, endosymbiont-host RNA-RNA interactions may drive maintenance of the symbiosis across fluctuating ecological conditions.

Expanded host and geographic range of tadpole associations with the Severe Perkinsea Infection group.

Severe Perkinsea infection is an emerging disease of amphibians, specifically tadpoles. Disease presentation correlates with liver infections of a subclade of Perkinsea (Alveolata) protists, named Pathogenic Perkinsea Clade (PPC). Tadpole mortality events associated with PPC infections have been reported across North America, from Alaska to Florida. Here, we investigate the geographic and host range of PPC associations in seemingly healthy tadpoles sampled from Panama, a biogeographic provenance critically affected by amphibian decline. To complement this work, we also investigate a mortality event among Hyla arborea tadpoles in captive-bred UK specimens. PPC SSU rDNA was detected in 10 of 81 Panama tadpoles tested, and H. arborea tadpoles from the UK. Phylogenies of the Perkinsea SSU rDNA sequences demonstrate they are highly similar to PPC sequences sampled from mortality events in the USA, and phylogenetic analysis of tadpole mitochondrial SSU rDNA demonstrates, for the first time, PPC associations in diverse hylids. These data provide further understanding of the biogeography and host range of this putative pathogenic group, factors likely to be important for conservation planning.

Characterization of the RNA-interference pathway as a tool for reverse genetic analysis in the nascent phototrophic endosymbiosis, Paramecium bursaria.

Endosymbiosis was fundamental for the evolution of eukaryotic complexity. Endosymbiotic interactions can be dissected through forward- and reverse-genetic experiments, such as RNA-interference (RNAi). However, distinguishing small (s)RNA pathways in a eukaryote-eukaryote endosymbiotic interaction is challenging. Here, we investigate the repertoire of RNAi pathway protein-encoding genes in the model nascent endosymbiotic system, Paramecium bursaria-Chlorella spp. Using comparative genomics and transcriptomics supported by phylogenetics, we identify essential proteome components of the small interfering (si)RNA, scan (scn)RNA and internal eliminated sequence (ies)RNA pathways. Our analyses reveal that copies of these components have been retained throughout successive whole genome duplication (WGD) events in the Paramecium clade. We validate feeding-induced siRNA-based RNAi in P. bursaria via knock-down of the splicing factor, u2af1, which we show to be crucial to host growth. Finally, using simultaneous knock-down 'paradox' controls to rescue the effect of u2af1 knock-down, we demonstrate that feeding-induced RNAi in P. bursaria is dependent upon a core pathway of host-encoded Dcr1, Piwi and Pds1 components. Our experiments confirm the presence of a functional, host-derived RNAi pathway in P. bursaria that generates 23-nt siRNA, validating the use of the P. bursaria-Chlorella spp. system to investigate the genetic basis of a nascent endosymbiosis.

A novel duplex qPCR assay for stepwise detection of multiple Perkinsea protistan infections of amphibian tissues.

Alveolate protists within the phylum Perkinsea have been found to infect amphibians across a broad taxonomic and geographic range. Phylogenetic analysis has suggested the existence of two clades of amphibian Perkinsea: a putatively non-pathogenic clade linked to asymptomatic infections of tadpoles in Africa, Europe and South America, and a putatively pathogenic clade linked to disease and mass mortality events of tadpoles in North America. Here, we describe the development of a duplex TaqMan qPCR assay to detect and discriminate between rDNA sequences from both clades of Perkinsea in amphibian tissues. The assay uses a single primer pair to target an 18S small subunit (SSU) ribosomal RNA (rRNA) gene region shared between the two clades, and two dual-labelled probes to target a region within this fragment that is diagnostic for each clade. This assay enables rapid screening for each of the two Perkinsea groups, allowing for detection, primarily of the phylogenetic group associated with disease outbreaks, and secondarily for the phylogenetic group with no current disease relationship identified. Incorporation of our novel qPCR assay into the routine surveillance of amphibian populations will allow for the assessment of the incidence of each protist clade, thereby providing an improved understanding of Perkinsea infection pervasiveness and a method to underpin future conservation planning.

A functional bacteria-derived restriction modification system in the mitochondrion of a heterotrophic protist.

The overarching trend in mitochondrial genome evolution is functional streamlining coupled with gene loss. Therefore, gene acquisition by mitochondria is considered to be exceedingly rare. Selfish elements in the form of self-splicing introns occur in many organellar genomes, but the wider diversity of selfish elements, and how they persist in the DNA of organelles, has not been explored. In the mitochondrial genome of a marine heterotrophic katablepharid protist, we identify a functional type II restriction modification (RM) system originating from a horizontal gene transfer (HGT) event involving bacteria related to flavobacteria. This RM system consists of an HpaII-like endonuclease and a cognate cytosine methyltransferase (CM). We demonstrate that these proteins are functional by heterologous expression in both bacterial and eukaryotic cells. These results suggest that a mitochondrion-encoded RM system can function as a toxin-antitoxin selfish element, and that such elements could be co-opted by eukaryotic genomes to drive biased organellar inheritance.

DivIVA Controls Progeny Morphology and Diverse ParA Proteins Regulate Cell Division or Gliding Motility in Bdellovibrio bacteriovorus.

The predatory bacterium B. bacteriovorus grows and divides inside the periplasm of Gram-negative bacteria, forming a structure known as a bdelloplast. Cell division of predators inside the dead prey cell is not by binary fission but instead by synchronous division of a single elongated filamentous cell into odd or even numbers of progeny cells. Bdellovibrio replication and cell division processes are dependent on the finite level of nutrients available from inside the prey bacterium. The filamentous growth and division process of the predator maximizes the number of progeny produced by the finite nutrients in a way that binary fission could not. To learn more about such an unusual growth profile, we studied the role of DivIVA in the growing Bdellovibrio cell. This protein is well known for its link to polar cell growth and spore formation in Gram-positive bacteria, but little is known about its function in a predatory growth context. We show that DivIVA is expressed in the growing B. bacteriovorus cell and controls cell morphology during filamentous cell division, but not the number of progeny produced. Bacterial Two Hybrid (BTH) analysis shows DivIVA may interact with proteins that respond to metabolic indicators of amino-acid biosynthesis or changes in redox state. Such changes may be relevant signals to the predator, indicating the consumption of prey nutrients within the sealed bdelloplast environment. ParA, a chromosome segregation protein, also contributes to bacterial septation in many species. The B. bacteriovorus genome contains three ParA homologs; we identify a canonical ParAB pair required for predatory cell division and show a BTH interaction between a gene product encoded from the same operon as DivIVA with the canonical ParA. The remaining ParA proteins are both expressed in Bdellovibrio but are not required for predator cell division. Instead, one of these ParA proteins coordinates gliding motility, changing the frequency at which the cells reverse direction. Our work will prime further studies into how one bacterium can co-ordinate its cell division with the destruction of another bacterium that it dwells within.

Dual Predation by Bacteriophage and Bdellovibrio bacteriovorus Can Eradicate Escherichia coli Prey in Situations where Single Predation Cannot.

Bacteria are preyed upon by diverse microbial predators, including bacteriophage and predatory bacteria, such as Bdellovibrio bacteriovorus While bacteriophage are used as antimicrobial therapies in Eastern Europe and are being applied for compassionate use in the United States, predatory bacteria are only just beginning to reveal their potential therapeutic uses. However, predation by either predator type can falter due to different adaptations arising in the prey bacteria. When testing poultry farm wastewater for novel Bdellovibrio isolates on Escherichia coli prey lawns, individual composite plaques were isolated containing both an RTP (rosette-tailed-phage)-like-phage and a B. bacteriovorus strain and showing central prey lysis and halos of extra lysis. Combining the purified phage with a lab strain of B. bacteriovorus HD100 recapitulated haloed plaques and increased killing of the E. coli prey in liquid culture, showing an effective side-by-side action of these predators compared to their actions alone. Using approximate Bayesian computation to select the best fitting from a variety of different mathematical models demonstrated that the experimental data could be explained only by assuming the existence of three prey phenotypes: (i) sensitive to both predators, (ii) genetically resistant to phage only, and (iii) plastic resistant to B. bacteriovorus only. Although each predator reduces prey availability for the other, high phage numbers did not abolish B. bacteriovorus predation, so both predators are competent to coexist and are causing different selective pressures on the bacterial surface while, in tandem, controlling prey bacterial numbers efficiently. This suggests that combinatorial predator therapy could overcome problems of phage resistance.IMPORTANCE With increasing levels of antibiotic resistance, the development of alternative antibacterial therapies is urgently needed. Two potential alternatives are bacteriophage and predatory bacteria. Bacteriophage therapy has been used, but prey/host specificity and the rapid acquisition of bacterial resistance to bacteriophage are practical considerations. Predatory bacteria are of interest due to their broad Gram-negative bacterial prey range and the lack of simple resistance mechanisms. Here, a bacteriophage and a strain of Bdellovibrio bacteriovorus, preyed side by side on a population of E. coli, causing a significantly greater decrease in prey numbers than either alone. Such combinatorial predator therapy may have greater potential than individual predators since prey surface changes selected for by each predator do not protect prey against the other predator.

Controlled sampling of ribosomally active protistan diversity in sediment-surface layers identifies putative players in the marine carbon sink.

Marine sediments are one of the largest carbon reservoir on Earth, yet the microbial communities, especially the eukaryotes, that drive these ecosystems are poorly characterised. Here, we report implementation of a sampling system that enables injection of reagents into sediments at depth, allowing for preservation of RNA in situ. Using the RNA templates recovered, we investigate the 'ribosomally active' eukaryotic diversity present in sediments close to the water/sediment interface. We demonstrate that in situ preservation leads to recovery of a significantly altered community profile. Using SSU rRNA amplicon sequencing, we investigated the community structure in these environments, demonstrating a wide diversity and high relative abundance of stramenopiles and alveolates, specifically: Bacillariophyta (diatoms), labyrinthulomycetes and ciliates. The identification of abundant diatom rRNA molecules is consistent with microscopy-based studies, but demonstrates that these algae can also be exported to the sediment as active cells as opposed to dead forms. We also observe many groups that include, or branch close to, osmotrophic-saprotrophic protists (e.g. labyrinthulomycetes and Pseudofungi), microbes likely to be important for detrital decomposition. The sequence data also included a diversity of abundant amplicon-types that branch close to the Fonticula slime moulds. Taken together, our data identifies additional roles for eukaryotic microbes in the marine carbon cycle; where putative osmotrophic-saprotrophic protists represent a significant active microbial-constituent of the upper sediment layer.

Unexpected mitochondrial genome diversity revealed by targeted single-cell genomics of heterotrophic flagellated protists.

Most eukaryotic microbial diversity is uncultivated, under-studied and lacks nuclear genome data. Mitochondrial genome sampling is more comprehensive, but many phylogenetically important groups remain unsampled. Here, using a single-cell sorting approach combining tubulin-specific labelling with photopigment exclusion, we sorted flagellated heterotrophic unicellular eukaryotes from Pacific Ocean samples. We recovered 206 single amplified genomes, predominantly from underrepresented branches on the tree of life. Seventy single amplified genomes contained unique mitochondrial contigs, including 21 complete or near-complete mitochondrial genomes from formerly under-sampled phylogenetic branches, including telonemids, katablepharids, cercozoans and marine stramenopiles, effectively doubling the number of available samples of heterotrophic flagellate mitochondrial genomes. Collectively, these data identify a dynamic history of mitochondrial genome evolution including intron gain and loss, extensive patterns of genetic code variation and complex patterns of gene loss. Surprisingly, we found that stramenopile mitochondrial content is highly plastic, resembling patterns of variation previously observed only in plants.

A single-cell genome reveals diplonemid-like ancestry of kinetoplastid mitochondrial gene structure.

Euglenozoa comprises euglenids, kinetoplastids, and diplonemids, with each group exhibiting different and highly unusual mitochondrial genome organizations. Although they are sister groups, kinetoplastids and diplonemids have very distinct mitochondrial genome architectures, requiring widespread insertion/deletion RNA editing and extensive trans-splicing, respectively, in order to generate functional transcripts. The evolutionary history by which these differing processes arose remains unclear. Using single-cell genomics, followed by small sub unit ribosomal DNA and multigene phylogenies, we identified an isolated marine cell that branches on phylogenetic trees as a sister to known kinetoplastids. Analysis of single-cell amplified genomic material identified multiple mitochondrial genome contigs. These revealed a gene architecture resembling that of diplonemid mitochondria, with small fragments of genes encoded out of order and or on different contigs, indicating that these genes require extensive trans-splicing. Conversely, no requirement for kinetoplastid-like insertion/deletion RNA-editing was detected. Additionally, while we identified some proteins so far only found in kinetoplastids, we could not unequivocally identify mitochondrial RNA editing proteins. These data invite the hypothesis that extensive genome fragmentation and trans-splicing were the ancestral states for the kinetoplastid-diplonemid clade but were lost during the kinetoplastid radiation. This study demonstrates that single-cell approaches can successfully retrieve lineages that represent important new branches on the tree of life, and thus can illuminate major evolutionary and functional transitions in eukaryotes. This article is part of a discussion meeting issue 'Single cell ecology'.

Environment-dependent fitness gains can be driven by horizontal gene transfer of transporter-encoding genes.

Many microbes acquire metabolites in a "feeding" process where complex polymers are broken down in the environment to their subunits. The subsequent uptake of soluble metabolites by a cell, sometimes called osmotrophy, is facilitated by transporter proteins. As such, the diversification of osmotrophic microorganisms is closely tied to the diversification of transporter functions. Horizontal gene transfer (HGT) has been suggested to produce genetic variation that can lead to adaptation, allowing lineages to acquire traits and expand niche ranges. Transporter genes often encode single-gene phenotypes and tend to have low protein-protein interaction complexity and, as such, are potential candidates for HGT. Here we test the idea that HGT has underpinned the expansion of metabolic potential and substrate utilization via transfer of transporter-encoding genes. Using phylogenomics, we identify seven cases of transporter-gene HGT between fungal phyla, and investigate compatibility, localization, function, and fitness consequences when these genes are expressed in Saccharomyces cerevisiae Using this approach, we demonstrate that the transporters identified can alter how fungi utilize a range of metabolites, including peptides, polyols, and sugars. We then show, for one model gene, that transporter gene acquisition by HGT can significantly alter the fitness landscape of S. cerevisiae We therefore provide evidence that transporter HGT occurs between fungi, alters how fungi can acquire metabolites, and can drive gain in fitness. We propose a "transporter-gene acquisition ratchet," where transporter repertoires are continually augmented by duplication, HGT, and differential loss, collectively acting to overwrite, fine-tune, and diversify the complement of transporters present in a genome.

Evolutionary diversification of the RomR protein of the invasive deltaproteobacterium, Bdellovibrio bacteriovorus.

Bdellovibrio bacteriovorus is a predatory deltaproteobacterium that encounters individual Gram-negative prey bacteria with gliding or swimming motility, and then is able to invade such prey cells via type IVa pilus-dependent mechanisms. Movement control (pili or gliding) in other deltaproteobacteria, such as the pack hunting Myxococcus xanthus, uses a response regulator protein, RomRMx (which dynamically relocalises between the cell poles) and a GTPase, MglAMx, previously postulated as an interface between the FrzMx chemosensory system and gliding or pilus-motility apparatus, to produce regulated bidirectional motility. In contrast, B. bacteriovorus predation is a more singular encounter between a lone predator and prey; contact is always via the piliated, non-flagellar pole of the predator, involving MglABd, but no Frz system. In this new study, tracking fluorescent RomRBd microscopically during predatory growth shows that it does not dynamically relocalise, in contrast to the M. xanthus protein; instead having possible roles in growth events. Furthermore, transcriptional start analysis, site-directed mutagenesis and bacterial two-hybrid interaction studies, indicate an evolutionary loss of RomRBd activation (via receiver domain phosphorylation) in this lone hunting bacterium, demonstrating divergence from its bipolar role in motility in pack-hunting M. xanthus and further evolution that may differentiate lone from pack predators.

Specialized proteomic responses and an ancient photoprotection mechanism sustain marine green algal growth during phosphate limitation.

Marine algae perform approximately half of global carbon fixation, but their growth is often limited by the availability of phosphate or other nutrients1,2. As oceans warm, the area of phosphate-limited surface waters is predicted to increase, resulting in ocean desertification3,4. Understanding the responses of key eukaryotic phytoplankton to nutrient limitation is therefore critical5,6. We used advanced photo-bioreactors to investigate how the widespread marine green alga Micromonas commoda grows under transitions from replete nutrients to chronic phosphate limitation and subsequent relief, analysing photosystem changes and broad cellular responses using proteomics, transcriptomics and biophysical measurements. We find that physiological and protein expression responses previously attributed to stress are critical to supporting stable exponential growth when phosphate is limiting. Unexpectedly, the abundance of most proteins involved in light harvesting does not change, but an ancient light-harvesting-related protein, LHCSR, is induced and dissipates damaging excess absorbed light as heat throughout phosphate limitation. Concurrently, a suite of uncharacterized proteins with narrow phylogenetic distributions increase multifold. Notably, of the proteins that exhibit significant changes, 70% are not differentially expressed at the mRNA transcript level, highlighting the importance of post-transcriptional processes in microbial eukaryotes. Nevertheless, transcript-protein pairs with concordant changes were identified that will enable more robust interpretation of eukaryotic phytoplankton responses in the field from metatranscriptomic studies. Our results show that P-limited Micromonas responds quickly to a fresh pulse of phosphate by rapidly increasing replication, and that the protein network associated with this ability is composed of both conserved and phylogenetically recent proteome systems that promote dynamic phosphate homeostasis. That an ancient mechanism for mitigating light stress is central to sustaining growth during extended phosphate limitation highlights the possibility of interactive effects arising from combined stressors under ocean change, which could reduce the efficacy of algal strategies for optimizing marine photosynthesis.